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To successfully carry out a wind energy conversion system, it is necessary to 
simultaneously control the rotor and the grid side. This paper proposes a 
doubly-fed induction generator's predictive power control. While the powers 
are controlled indirectly through currents, the latter is controlled using the 
deadbeat command. Based on discrete-time, the control suggests at each 
sample period the required voltages to the the-back-to-back converter to 
reach the desired setpoints, control the powers, and the DC link voltages. For 
these reasons, a presentation of the system is given first, then a description 
of the predictive control, followed by applying this strategy on the rotor side 
control and grid side control. Finally, a random wind profile was applied to 
analyze the system's performance with a unitary power factor. The 
simulation results are presented in the MATLAB/Simulink environment 
using a 1.5 kW DFIG. The results obtained by applying a random wind 
profile have well fulfilled the objectives of the control and the system 
robustness is approved by the excellent tracking allowing the machine's 


internal parameters variation. By comparing the quality and the tracking 
reference of the proposed control method to other control methods, the 
deadbeat controller was very promising. 
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1. INTRODUCTION 

The energy transition is defined as the subtraction of fossil energy to clean energy. To limit the 
climate change effect, greenhouse gas (GHG) emissions need to be reduced in the energy production sector 
[1]. According to the international renewable energy agency (IRENA), the promoter's actions and measures 
of renewable energy and energy efficiency could reach 90% of the required carbon reductions by 2050. For 
that purpose, control information and smart technologies need to be more developed [2]. 

To deal with this situation, Morocco has adopted an ambitious national strategy to reduce GHG 
emissions and secure its electricity supply through the large-scale development of renewable energies and 
energy efficiency. For that purpose, significant reform of the legislative and regulatory framework has been 
initiated. This profound change in the energy sector towards renewable energy must be accompanied by a 
change in the practices and technologies of generation, transport, and energy use. 

In fact, Morocco had 3000 MW installed renewable energy (RE) capacity (1770 MW of hydro, 1015 
MW of wind, and 180 MW of solar) according to the Moroccan agency for renewable energies (MASEN) in 
2017. Today, Morocco faces a new challenge: the Kingdom has decided to increase clean energy in the 
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electricity mix from 33% in 2017 to 42% by 2021 [2]. Within this framework, the Moroccan program for the 
development of wind energy falls within the framework of public/private partnerships or private projects of 
new wind farms, which will bring the installed electrical power of wind power to 2000 MW by the end of 
2021. 

In this context, this work aims to elaborate an optimal design of a wind power system. By looking at 
the literature, Kumar & al. proposed a review of the different wind turbine technologies and three kinds of 
wind systems were presented: Permanent magnet synchronous generator (PMSG), squirrel cage induction 
machine (SCIG), and doubly fed induction generator (DFIG) [3]. The most common promoters of the wind 
turbine industry, such as siemens Gamesa, Suzlon, and Enercon, commonly use DFIG with installed powers 
ranging from 1.5MW to 6MW to produce wind energy conversion systems (WECS). DFIG has the most 
significant average power capacity [4], thanks to its ability to generate power from both stator and rotor. This 
advantage came from its capacity to operate into two modes: sub synchronism and hyper synchronism, 
depending on wind speed [5] . 

For practical, stable, and optimal energy production, active and reactive powers, currents, and 
voltages need to be controlled. Among the control strategies, vector control is used widely. To simplify the 
control, field-oriented control (FOC) was used to make DFIG similar to direct current (DC) via the 
proportional integrator (PI) regulator. Therefore, PI controllers have some drawbacks, such its limited 
performance under external disturbance [6]. To deal with this issue, other alternatives have been proposed. In 
[7] authors proposed adaptative backstepping control as an alternative to FOC. This control uses the 
Lyapunov function to ensure better system stability [8], [9]. Nevertheless, on the other hand, determining the 
Lyapunov function is complicated because it relies on experience and luck. Arrousi et al. [10], the authors 
suggested a detailed direct power control (DPC) study based on hysteresis comparators. DPC's main 
advantage is its independence from system parameters, and decoupling between system variables is no longer 
needed [11]. In [12] authors suggested sliding mode control. This technique deals very well with 
uncertainties and parameter variations to make the system more robust [13]. However, this latter is the origin 
of chattering phenomena [14]. 

On the other hand, predictive control is an effective control strategy for machines and converters as 
an alternative to FOC [15]. It is based on the system model, which can predict the behavior of controlled 
variables. It suggests an optimal control for each sampling interval. It is exceptionally efficient since it can 
manage the constraints, system states, and output variables [16]. However, some researchers in predictive 
control focused either on the rotor side control (RSC) [17] or the grid side control (GSC) aside [18]. 
Although to successfully carry out the system, it is necessary to control rotor and Grid Side simultaneously. 
In the interest of filling this gap, the deadbeat predictive controller is applied to control the RSC and the GSC 
of the DFIG. The main purpose of this control is to ensure a maximum of active power (tracking the shape of 
wind) with a unit power factor sinusoidal currents and voltages; with a normalized frequency and a low total 
harmonic distortion (THD) to be able to inject the produced electricity into the grid. This aim is supposed to 
be fully filed by ensuring the closest measured variables to the predictive ones that need to be applied to the 
back to back converters. 

In order to apply the necessary voltages, rotor currents need to be controlled through the predictive 
regulator. The rotor side converter provides the voltages required to control the power exchanged between 
the stator and the grid. The grid side converter controls the bi-directional power exchange between the grid 
and the rotor by adjusting the Dc link and grid currents. 

This paper aims to develop the predictive current control method for WECS. This method provides 
reference tracking of current and carries out power tracking by eliminating errors in each sampling period. It 
provides a faster dynamic response like conventional model predictive control methods. Besides, evaluating 
the switching state or voltage level, calculating the cost function, and selecting the weighting factor is not 
required. The developed control is designed to ensure optimal operation and improve the quality of energy 
produced by WECS through the combination of control between RSC and GSC and acting on the system 
performance and response time. 

The introduction is presented in the first section of this article. The mathematical modelling of the 
wind energy conversion system is discussed in section 2. Section 3 introduces the deadbeat method and 
provides a detailed mathematical demonstration of the proposed control applied to RSC and GSC. The 
simulations results are presented and discussed in section 5; a performance’s comparison is also provided 
with previous works. Finally, a conclusion is given in the last section. 


2. MATHEMATICAL MODELLING OF THE SYSTEM 
The WECS system is composed of a turbine linked to DFIG via a gearbox to adjust the speed. The 
stator is connected directly to the grid. Therefore, the rotor is linked to the grid by the back-to-back 
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converters. This latter has the task of applying the control to the system. WECS system consists of producing 
electrical energy from the wind. The power of the wind (Pv) produced by a turbine depends on three 
parameters: air density (p), wind speed (V), and the area swept by the blades. The power recoverable by a 
wind turbine is proportional to the surface swept by its rotor (S) and the cube of the wind speed [19], [20]. 


p= ose (1) 
2 
The DFIG stator is linked directly to the grid, while its rotor is linked via a back-to-back converter. Based on 
the DFIG space vector model, the mathematical model of the DFIG is expressed in the (d-q) reference frame 
as follows [21], [22] . Where, R, L refers respectively to Resistance and Inductance (machine's parameter). 
The indicators s and r refer to stator and rotor. The voltages in the referential (d, q): 


Vsa = Rs: Isa + 284 = sq. Os (2) 
Vaq = Rs: leq + 28 + Psa. Ws 6) 
Ved = Rr lng + 24 prg Or (4) 
Vea = Ryo Irq + 22 pra Or (5) 


The fluxes in the referential (d, q): 


Psa = Ls Isa + Lm. Ira (6) 
Psa = Ls- Isq + Lm: Irq (7) 
Ora = Ly. Ira + Lm. Isa (8) 
Org = Lr: Irq + Em: Isq (9) 


The active and reactive powers: 


Py = Vsa: Iga + Voq: Iq (10) 
Qs = Veq: Isa — Vsa- Isq (11) 
P, = Vra: Ira + Vra: Irq (12) 
Qr = Vra: Ira — Vra: Irq (13) 


The power converter is an essential component for controlling the WECS. It is composed of two converters 
(AC/DC/AC) linked by a DC bus. Its function is to regulate the DFIG's powers and currents, then provide the 
optimal active power and unit power factor to the grid [23], [24], [25]. 


3. DEADBEAT STRATEGY 

For a specific control variable, predictive control uses the system model to anticipate future process 
behavior then achieve optimal action based on predefined optimization criteria. This strategy identifies the 
appropriate control variables that cancel the error between the control variable and the reference input [18]. 
In general, the continuous representation of a system is presented as [26], [27]: 


x(t) = Ax(0) + Bu(t) + Gw(t) 
f y(t) = Cx(0) ay 


x (t) represents the stator vector to be controlled (e.g. voltage, flux, speed, current, torque, power...) 
u (t), y(t) represents input and output vectors, respectively (e.g. voltage, current, torque.) 
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w (t) perturbation 
A, B, G are n*n matrices and C Identity matrix 

The representation of the system in the discrete-time is necessary to know the behavior of the 
system. It is written as [28]: 


Xk+1 = Ax(k) + Bu(k) + Gw(k) (15) 
This discretization is obtained by making the following transformation: 

A=eT=1+AT,B= H e^TB dt = BT, G e eATG dt = GT (16) 
To obtain the deadbeat control, the reference is set equal to the next sampling instant: 

X k+1 = Xref (17) 


The (18) is obtained by injecting (17) in (15): 


i(k) = B~. X ref — A. B74. X(k) — B-t. G. wa(k) (18) 
The deadbeat control method operates under a high proportional controller gain. On the other hand, resulted 
power must track the reference by ensuring a steady-state null error. The input corrector can be written as 
(19): 

ū(k) = Fa — X) (19) 


Where F is the gain matrix. 
The requisite input is determined using the following formula: 


U(k) = B-t. AHA. Xref — X(k)] — B74. G. Wa(k) (20) 


3.1. Rotor side control 

By adjusting the rotor currents (Ird and Irq) through the predictive regulator, the rotor side converter 
(RSC) provides the voltages required (Vr) to control both Ps and Qs [29]. Vector control is applied to ensure 
decoupling between direct and quadrature components [30]. Therefore, the stator flux vector is aligned with 
the d axe. As results, the flux and voltage equations become: 


®sq =0, Psa = Qs Q1) 
Vsa = 0,Vsq = sd: Ws (22) 
In (23) and (24) are obtained by injecting (8) and (9) into (4) and (5) respectively. 


dira 


la (23) 


Ved = Ry- La = Ly: Ory. leq +0 


dl; 
Vra = Re: Irq — Lr- Or- Ira + 0.2 + Ly. Or: Isa (24) 


Using the previous equation, a state-space representation is given by the: 


dira —Rr wr.LM 

dt L o.Lr Vra o.Ly 

=f ESE f k ii ea Veal |=» Torla a =i lal (25) 
o.Lr. 


dt 


Lo” 
Ls-Lr 


with: o = 1 — 


By applying (14) on the system parameters, the rotor currents can be discretized as: 
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R,.T @y.T T @y.Ly.T 
TA 1- o.Ly o g _ [olr 0 0 o.Ly 
Ar rT j Ret B; T |r T |-or.LmT 0 (26) 
o o.Ly o.Ly o.Ly 
The resulting rotor voltages equations can be expressed as: 
Ly 
Vea(k) = E (Tra (k + 1) — Ip (9)+R y. Tra (k) — Lr- Or- Irak) — Lu: @y-I5q(K) (27) 
Ly 
Veq(k) = E (Inq ( + 1) = Irq K) ) + Rr- Ira 00 + Lr- @p. Ira (k) + Lu: 051545) (28) 


It can be seen from (27) and (28) that this control method handles the machine variables, including the stator 
currents. 


Z2 M 
Ps = =E. E. Ve- Irq (29) 
3 s 
Qs = —3.V; a) (30) 


Considering the assumption of (17), the current references (Ird_ref and Irq_ref) are set equal to the currents at 
the next sampling time Ird(k+1) and Irq(k+1) respectively. The Irq_ref and Ird_ref are expressed as: 


2 Pref.Ls 


lrq ref =r Vile (31) 
2 Qref-Ls s 
lra ref = EE y (32) 


Vs.Lm Lm 


3.2. Grid side control 

By controlling the grid currents Igd and Igq through the predictive regulator, the grid side converter 
(GSC) provides the voltages required (Vf) allowing the control of the active power bi-directionally 
exchanged between the rotor and the grid. The power factor and the DC bus voltage are controlled as well. 
By orienting the grid voltage (Vg) along the q axis, Vgd is set equal to zero Vgd=0 & Veq=Vg. The voltages 
exchanged between the Grid supply and Grid side converter in the continuous dq frame become: 


Vig = = Ry. Iga + ipo a — Le. Ws. l2q (33) 
dlgg 
Vig = Veq + Rs. Izq + Le. ETE + Le. Ws. Iga (34) 


Using the previous equation, state-space representation is given by the: 


diga 
dt alil A el, 9 = Le] (35) 
dlgq | I ral vo] + 0 Té Veq 


dt 


By applying (14) on the system parameters, the filter current can be discretized as: 
bank + 1) = Aglan 09 + Bg. Vrca q) 09 + Gp. Vga, (19 (36) 


The following discrete-time parameters are obtained by applying the required transformations given in (16) 


eT ee 0 0 
A,=| ” -T 
=> al ale 


—05. T Le 
The resulting voltages equations can be expressed as: 


Novel deadbeat predictive control strategy for DFIG’s back to back power converter (Manale Bouderbala) 


144 O ISSN: 2088-8694 


Va Œ) = E (Iga(k + 1) = Iga (k)) + Re- Tga (1) — Le- 005: Iga (K) (37) 
Viq(k) = E (Igq(k + 1) — Igq(K)) + Rr-Igq(9 + Lp. 05. Iga (K) + Vog (k) (38) 


Considering the prediction variables Igd(k+1) and Igq(k+1) represent the grid reference currents Igd_ref and 
Igq_ref respectively. And based on the vectorial approach used and the grid power expressions, the reference 
grid currents are given by: 


a Pg ref 

Igq_ref = Vea (39) 
_ Qg ref 

lgd ref = Veq (40) 


The current Igd_ ref is obtained directly from Qg ref while the current Igq_ref is obtained by controlling the 
DC link voltage. Figure 1 shows the whole system with the control principle. 


C_ref - 


ator 


Figure 1. Deadbeat control applied to WECS 


4. RESULTS AND DISCUSSION 

In order to test and validate the control developed previously, some tests have been performed. First 
of all, a random wind profile with a zero reactive power setpoint is applied to the system. Figure 2 shows the 
wind profile used, which has an average value of 8m/s, while Figure 3 shows the electromagnetic torque that 
varies inversely with the wind profile. Figure 4 and Figure 5 show the active and reactive stator powers. 

It is clear that the decoupling of these two powers is assured. The wind profile has the same inverse 
form as the active power, which follows its reference with a 20W error and a response time of 0.6 s. The fact 
that the equipment is in generator mode explains the negative figure of active power. The reactive power 
follows its reference with a response time of 0.4s and a deviation of (5Var). Figure 6 illustrates the 
quadrature rotor current. It is shaped like the inverse of active power as shown in Figure 2. Ps changes 
linearly with Irq through the negative coefficient indicated in (31) therefore this result was expected. Figure 7 
illustrates the direct rotor current. Since Qref is set to zero and the remaining constant ratio s/Lm (32) is equal 
to 5 A, it has a constant value of about 5A. Rotor quadrature and direct current errors were found to be fewer 
than 0.025 A. This little error proves that they effectively reached their references. 

As the current is the power picture, Figure 8 depicts stator currents with an inverse form matching to 
the wind profile. While the profile fluctuates, the currents maintain a sinusoidal form with a period of 0.02 s, 
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providing a frequency equal to the grid frequency of 50 hz. The normalized value of THD must be less than 
5% to have acceptable quality power. In this investigation, the THD was less than 0.27% as shown in Figure 
9, indicating good power quality. Figure 10 depicts the rotor current profile, which fluctuates sinusoidally 
and proportionately with the electromagnetic torque. From t=6 s to t=8 s, the current form transitions: the 
generator switches from subsynchronism to hypersynchronism mode. Concerning VDC response, it can be 
seen in Figure 11 that Dc link voltages need 0.18s to reach the constant reference value of 600 V. An 
overshoot of 8.3% explains the VDC peak. 

On the other hand, the grid active and reactive powers are illustrated in Figure 12 and Figure 13, 
respectively. The Pg has the same shape as the Ps profile, and Qg profile has the same shape as the Qs 
profile. As for the Ig, Figure 14 shows that the currents grid has the same sinusoidal shape and the same 
period as Is. These results testified in the interest of controlling RSC and GSC simultaneously tracking the 
shape of the wind to inject the produced electricity into the grid. As shown in Figure 15, the voltage and 
current of one phase of the grid are sinusoidal with the same period of 0.02s in phase opposition, ensuring the 
unit power factor. 

To test and verify the control’s robustness and performance, the ability of Ps to track the reference 
needs to be checked. A step has been applied as a reference, while stator and rotor resistances have been 
changed from Resistance/2 to Resistance*2. As shown in Figures 16 and Figure 17 can be seen that the 
tracking is still ensured despite the slight increase of the oscillations. A comparison of results between the 
developed control and other recent studies is shown in the Table 1. It should be mentioned that they do not 
refer to the same conditions since it is very difficult to find several works done under the same conditions. 
Although the response time is slightly higher than the response time proposed by [8], the biggest advantage 
of the deadbeat control is the significant reduction of the tracking error, and the error was reduced by 65% 
compared to the BSC and 87% compared to the FOC. Comparing the THD of this study with [31] and [12], it 
is remarkable that the THD has been reduced in a very significant way for the four types of controls 
proposed. All this leads, the Deadbeat controller proposes the best produced energy in terms of performance, 
robustness and quality. 
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Table 1. Performance's comparison 
Ref Control method Response time (s) _ Error (%) _ THD (%) 


[8] Backstepping 0.3 1.17 - 
Vector control 0.8 3.3184 - 
[31] Direct Power Control - -- 5.17 
Vector control - - 5.95 
[12] Backstepping - - 2.36 
Sliding mode - - 3.99 
[14] Sliding mode - - 2.98 
Deadbeat controller 0.4 0.4 0.27 


5. CONCLUSION 

The predictive power regulation of a doubly-fed induction generator was examined. The deadbeat 
command was designed to control the powers systems through the currents. The deadbeat control 
recommended the appropriate voltages to the back-to-back converter at each sample period based on discrete- 
time to achieve the required set-points, control the powers and the DC link voltages. To guarantee decoupling 
between the two components, axes d, and q, field-oriented control was used. For these goals, the system was 
provided initially, followed by a discussion of the deadbeat predictive control approach, and finally, applying 
this strategy to rotor and grid side control. With a unitary power factor, a random wind profile was utilized to 
examine the system's performance. Using a 1.5 kW Doubly Fed Induction Generator, the simulation results 
are shown in the MATLAB/Simulink environment. 

With the use of a random wind profile, the results produced have met the control and system 
objectives by assuring good reference tracking; they have also met the normalization standards, yielding a 
THD of 0.27% and a frequency of 50 hertz. The excellent tracking that allows the machine's internal 
characteristics to vary attests to the control's durability. The deadbeat controller was found to be quite 
promising as compared to other control systems in terms of quality and tracking reference. Once these results 
in the MATLAB software are produced, as a perspective, the application of these commands on the test 
bench is conceivable. 
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